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ABSTRACT 

Simplified models are developed for radiative heating and cooling 

s 

and for ozone photochemistry in the region 22-60 km. The latter permit 
the inclusion of nitrogen and hydrogen reactions in addition to simple 
oxygen reactions. The simplicity of the scheme facilitates the use Ox 
a wide variety of cooling and reaction rates. We also consider deter- 
mination of temperature and composition as a j oint' process . It is 
shown that joint radiative-photochemical equilibrium is appropriate to 
the mean state of the atmosphere between 35 and 60 km. Equilibrium 
calculations are then used to show that hydrogen reactions are important 
for ozone and temperature distributions primarily above 40 km while 
nitrogen reactions are important primarily below 50 km. Comparisons 
with observed distributions of temperature and ozone suggest the need 

for water vapor mixing ratios of from 0.5-2. 5 x 10 and mixing ratios 

-8 -7 

of (W0 2 + NO) of from 2 x 10 -1.6 x 10 at the s'tratopause . At 35 km, 

_8 

a mixing ratio of (N0 2 + NO) of about 3 x 10 is indicated. The 
precise values depend on our choice of reaction and radiative cooling 
rate coefficients, and the simple formulation permits the reader to 
check the effect of new rates as they become available. 

The relaxation of perturbations from joint radiative-photochemical 
equilibrium is also investigated. In all cases the coupling between 
temperature dependent ozone photochemistry and radiation lead to a 
reduction of the thermal relaxation time from its purely radiative 
value. The. latter, which amounts to about 10 days, is reduced to 2-4 
days at heights of 31-35 km. This greatly enhances the dissipation or 
waves travelling through the stratosphere. 
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1. INTRODUCTION 

The equilibrium temperature in the stratosphere is determined 
by the approximate balance between' heating due to absorption of solar 
.ultraviolet energy by ozone and cooling due to infrared emission by the 
15p band of carbon dioxide. (A certain amount of cooling is also due to 
O^ 9.6p emission, but this cooling is considerably less than that due to 
the 15p band of CO^ [Murgatroyd and Goody 1958]). The rate of energy 
absorption and therefore the temperature depend on the ozone density. 

Production of ozone is by the reaction 

O + O z + M 0 3 + M U; 

which is strongly temperature dependent. This coupling between tem- 
perature and ozone density indicates that the relaxation time of a tem- 
perature perturbation should not be that due to carbon dioxide cooling 
alone but should be that due to carbon dioxide cooling modified by the 
effects of photochemistry. 

Lindzen and Goody (1965) have calculated the thermal relaxation 
time for carbon dioxide cooling coupled with photochemical effects 
for a pure oxygen atmosphere. Since these calculations were made, 
the rates for pure oxygen reactions have been questioned and alter- 
native rates have been proposed. In addition it now seems that 
reactions involving nitrogen and hydrogen compounds may be sig- 
nificant in determining the ozone distribution in the stratosphere. 
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(Hunt 1966, Leovy 1969a, Crutzen 1911 ) and such reactions should there- 
fore be included in the photochemical model. Hoot of the rates for 
reactions involving hydrogen compounds are not temperature dependent 
so there is a question as to the influence of these reactions on the 
coupling between prone density and temperature. It has been suggested 
(Leovy 1969a) that if the hydrogen reactions are dominant in determining 
the ozone density then the appropriate thermal relaxation time is just 

this is untrue since relation (1) remains important in all schemes. 

Cooling due to infrared emission by carbon dioxide and modified by 
photochemistry acts as damping on motions in the stratosphere, 
time scale for this damping, which may be represented as Newtonian 

cooling, is simply the thermal relaxation time scale. As Mchinson 

i- ' excited in the troposphere, could 

(1969) shows, internal Rossby waves, exerted in 

carry large amounts of energy up to the lower atmosphere unless this 
time is much less than ten days, ten days being the approximate time 
ecale for unmodified C0 2 15 m cooling. Thus, the effects of photo- 
chemistry on this time scale could be crucial. 


In this paper we shall investigate the effect of photochemistry on 
the cooling rate for the height range 22-61 km. A simplified photo- 
chemical model that is suitable for this height range and that in- 
cludes reactions involving nitrogen and hydrogen compounds is de- 
veloped in section 2. The temperature equation for this photo- 
chemical model is obtained in section 3. Vertical distributions of 
constituents and temperature for radiative-photochemical equilibrium 
are discussed in section 4. Linearized equations for the photochemistry 
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and temperature are found in section 5 as are the relevant time scales 
for photochemical and thermal relaxation. Conclusions arc presented 


in section 6. 

Many reaction rates are uncertain by as much as an orde. of 
magnitude so a wide range of values is used. Because mixing ratios 
for water vapor and nitrogen oxides are uncertain several distribu- 
tions are used. For the range of rates and ratios used, it is found 
that the effect of including photochemistry is to reduce the thermal 
relaxation time scale above 30 km from the value appropriate when 
only infrared emission by carbon dioxide is considered. Further- 
more, the minimum value for the thermal relaxation time scale 
always occurs at heights between 31 and 35 km and with values be- 


tween two and four days. 

* It is also found that equilibrium distributions of temperature and 
ozone density vary with the different rates and ratios assumed. At 
the stratopause relaxation time scales are short enough (^90 
minutes for ozone and a few days for temperature) that the equilibri- 
um temperature should be close to the observed temperature of ap- 
proximately 270°K (Supplemental Atmosphere [Tropical], Theon and 
Smith 1971). We attempt to match the observed temperature rather 
than the observed ozone densities since the former is better deter- 


mined than the latter between 35 and 61 km, 
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The neglect of nitrogen and hydrogen reactions leads to the predic- 
tion of excessive temperatures between 35 and 61 km. The prediction of 
observed temperatures between 50 and 61 km calls for water vapor mixing 
ratios of from 0.5-3 x 10 ^ depending on the reaction rates used. A 

O 

mixing ratio for total nitrogen oxides (NO + NO 2 ) of 2-3 x 10 produces 
correct temperatures from 35-50 km. All such estimates depend on our 
choice of reaction and cooling rate coefficients, but the dependence 
has been specified in a simple manner. 
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2. PHOTOCHEMICAL MODEL 
As many as 80 different reactions have been considered for an 
atmosphere consisting of nitrogen, hydrogen, oxygen, and their 
compounds (Hunt 1966, Shimazaki and.Lai-rd 1970, Crutzen 1971). 

a . 

For many problems the full set of reactions may not be needed and a 
much shorter set can be used for the photochemical model. It should 
be apparent that the simplified photochemical model developed in this 
section may not be the most suitable one for all problems. 

For the current problem we are interested primarily in the 
vertical distribution of ozone, the reactions which determine this 
distribution, and the temperature dependence of such reactions for the 
height range, 22-61 km. Many reactions are important in only limited 
height ranges. Reactions which are important only above 61 km or 
below 22 km are not retained. The list of reactions retained is in 

Table 1. 

For each constituent found in Table 1 it is possible to write a 
continuity relation 

dn. . 

— - + n. V * v = S. - L. (2) 

dt i l i 

where d/dt is the substantive derivative, -^ L + , ~v is the velocity 

vector and n S , and L. are the number density, source term and 
* i’ i i 

3 1 

loss term, respectively, of the constituent i(i = 0(0( F)), 0*(0( D}), 
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70 km 


n./» ■ (( 1 

i m 


•where 


B n is the molecular number 


density, so the continuity equation 


for n becomes 
m 


-^-n +n V * v - 0- 

dt m m 


Substituting eq (4) mto eq ( 2 ) yi elds 


dn./n 

__i — 22- = (S. - L-)/ n T 

At l i 1 


The ratio n./n^or a constituent i is referred to as-its mixing ratio. 

xhe photochemical model now consists of 13 continuity equations 
of the form (5). There are several relations among the constituents 

some of these relations have 

which will further simplify the model. 

i mYitsch 1968, Leovy 1969a, and Crutzen 
been pointed out previously (Dutsch 


1971). 

Since below 61 km less than 1% of the H.O molecules are is- 
sociated in a day by reactions (3, and (8), the mixing ratio of water 
vapor at a given height may be considered constant. Also the^ 

adjustment time for equilibrium of O* is t(0*> = <Vm> < 10 

. ' i, ( 4' sec. 

and. that for H is.t(H) = ( n Q m 10 


Both of these 
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constituents may be considered in equilibrium with concentrations 


q 3b n 0 3 

n 0* # k 7 n 

7 m 

and 

^ n QH n 0 k 9 

“n ^ n 0 2 n m k 10 + n 0 3 k 27 


( 6 ) 


(7) 


As Leovy (1969a) indicates, eq (7) means that when reaction (27) is 
small compared to reaction (10), the atomic hydrogen produced in 
reaction (9) is destroyed so rapidly by reaction (10) that the two 
reactions may be replaced by the single reaction 

OH + 0 — > It0 2 


with the reaction rate, k 9 . This approximation is. valid below 40 km 
but about 5% of the hydrogen destruction at 40 km is by reaction (27) 
and above this height reaction (27) becomes increasingly important to 
the destruction of hydrogen and ozone, accounting for 20% at 60 km. 
Thus, if reaction rates are significantly alterred, reaction (27) could 
become a major one. 



When the reactions involving 0 and 0^ are considered, it is 
apparent that the dominant reactions involve interconversion between 

0 and 0^. The ratio 


rv 


q 3a + q 


3b 


n _ n k,. 
0 2 m 5 


has an adjustment time, tCRj) - (q 3a + q 3b + n 0 “mV < 10 S£C> 

which is short compared to other time scales of interest so eq (9) 
remains valid. Similarly the reactions involving .interconversa.cn 
between 110 2 - and OH are dominant among those involving these two con 


stituents . 


Below 40 km, where (8) may be used 



11 


R- 


n 


’HO, "o k 9 +n 0 3 k 12 

1 r±t. — 


OH 


~ n O k ll + n NO k 15 


( 10 ) 


which has an adjustment time of t (R^) “ ( n 0 k 9 + n o 3 k 12 + n 0 k ll + 

n k < 100 sec. Above 40 km 

NO 15' 


n 


OH 


zy 


“U0 2 n 0 k ll + n H I ’Oj k 27 

"A 


(ID 


and 


"HO, ~ 


n H n 0 2 n m k 10 + n 0H n P 3 k 12 
n O k ll 


(12) 


which have adjustment times t(0H) “ (n^) < 15 sec and 

(jiok^i)" 1 < 15 sec. Therefore, when eq (7) is used in eq (12) the 

resulting ratio is . 


n 


Rn 


ho 2 k g 

n 0H ~ k ll 


"° 3 k 21 .! 

C1 + n n n kT 5 
0 2 m.10 


(13) 


+ k 


12 


R l k ll 


The term — in (13) becomes negligibly small above 45 km. 

Vu 

Finally, by similar arguments, one gets a ratio, 


n. 


R, 


NO 


n 


rv 


NO. 


n 0 k 16 + . q N0 2 
n 0 3 k 17 


(14) 
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+ n k ) _1 ( 200 sec 

which lias an adjustment time t(Rj) = ( n 0 k 16 + V0 2 ° 3 17 

UNO HNO , and NO are approximate.! 

The three constituents, HN0 3> 2 3 

in equilibrium with values of 

• VmNo/oH (15) 


n 


HNO, 


k 2i n O 


n 


HNO, 


k ?.n n m n OH n NO 


k 22 n O 


( 16 ) 


and 


n 


NO, 


n 


k Zi n o n HNO, + k 18 n N0 2 "0 3 


k 23 n NO 


(17) 


• t(HNO ) = (k 7 ,n f 1 < 2 * 10 sec ' 

Their adjustment times are. T(HN0 3 ) l 21 Q ^ 3 


3 -andT<N0 3 >=( k 2 3 n N o r <2 ' 10 ^ 

22 the use of equilibrium value 


t(HNO ; ) = (Vo) (2 ' 10 560 


These times may 


be sufficiently long for 
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to be questioned. However, eqs (15), (16), and (17) permit an estimate 

of the densities of HN^, HNO ? and NOj relative to the densities of NO, 

NO HO and HO,. These estimates indicate that the constituents HNO^, 
2 ' 2 

HNO , and NO are of minor importance in determining the vertical 
2 3 

distribution of O . Therefore eqs (15), (16), and (17) are valid approxi- 
mations for the present problem. 

The total odd nitrogen density is the sum of the densities of NO, 

NO NO , HNO_, and HNO . Since the list of reactions in Table 1 
2 3 Z 3 

does not contain sources or losses for the total odd nitrogen density, 
the continuity equation is 


_d 

dt 


t n NO + n NO z + n N0 3 + n HN° 3 ‘ n HNO z J 
n 

m 


(IS) 


or, using eqs (15), (16), and (17) in eq (18), one gets 


jd_ 

dt 


»NO + n NO z ] 
n 

m 



(19) 


Two sources for total odd nitrogen in the region 22-61 km are con- 
sidered to be downward diffusion of NO from above 70 km where it 
is formed and the upward diffusion of from the troposphere, 

which is then destroyed by the reaction 
N z O +. O* 2NO 


( 20 ) 
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-11 3 -1 

with a rate k = 9x10 cm sec (Donovan and Husain 19?0). Now' the 
• . . - 7 

mixing ratio of N^O is 2. 5x10 at the tropopause and apparently de- 

, s 

creases with height (Schutz et al. 1970). If we estimate the mixing 
ratio of NO + NO^ as 10 - 10 (Junge 1963, Barth 1966, Pearce > 

1969), then the time scale for an increase in odd nitrogen density due 
to the two sources described is estimated to be several weeks. The 
possible losses to the odd nitrogen density are downward diffusion of 
the constituents or conversion to N^O^. £°ll°wed by downward diffusion 
of that constituent. 


Rather than try to include the source and loss terms for the odd 


niti ogen concentration, we have decided to use eq (18) or, equivalently, 
to set 


( "NO + "no. 


)/n = K 


m 


( 21 ) 


where K is a mixing ratio which is constant in time but a possible 
function of height. This highly simplified treatment of the nitrogen 
constituents may be justified on the grounds that we are not seeking 
a detailed time and height distribution for e ach nitrogen constituent 
but are trying to evaluate the importance of the nitrogen reactions 
t o the distributi on of ozone. 

The constituent H^O^, is produced in reaction (24) and destroyed 
in reactions (25) and (26). The net result is the loss of IIO back 
into H^O. This conversion also occurs in reaction (14). If reactions 
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(24) , (25), and (26) are neglected the result is that the sum, n OH + n HO * 

2 

calculated is too large. However, the change in the sum, n Qj.j + 
occurring when these reactions are ignored, is significant only below 
35 km. It will be shown in sections 4 and 5 that the reactions involving 
hydrogen are significant only above 40 km. Therefore, reactions (24), 

(25) , and (26) are to be neglected. 

The constituents O*, H, HNO^> HNO^i and NO^ are considered to 
be in equilibrium (eqs (6), (7), (15), (16), and (17) and the water vapor 
mixing ratio is considered constant in time while H-^O^ is neglected. 
Therefore, the continuity relations, eq (3), are needed for only the 
six constituents, O, O , NO, NO^ HO, and HO^ By using eqs (9), 

(10), and (14), the six continuity relations are reduced to three: 



( 22 ) 

(23) 

(24) 


with 
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C «* 2 n n /n 
<1 2 0 2 m 

f - “iVm 




B - (Rjkg + R 1 H 2 k n + k 12 


k,On z < AO km 


m 


n„ n k_ n 
0 m 10 ]_ 


B « (R x k 9 + R 1 R 2 k ll + R l k 9. (1 + n 0 k 27 } 


+ k.,)n m z > AO km 
12 ra 

F =■ (2 R 1 k l6 + k l8^ n m 
D » 2k 8 n H 2 O q 3b^ k 7 n m 
G “ 2q H 2 0 n H 2 0 /n m 
E - 2(k 13 + R 2 k i4 )n ra ' 


(25) 

(26) 
(27) 


(28) 

(29) 

(30) 

(31) 

(32) 


Using (13) in (28) we get 


B £ 2(R x k 9 + k 12 )n ra 


(33) 


„e formally retained (27) above 40 Urn, the me would also obtain 
) as an approximation. Thus (27) may be used at all heights with 
, small errors. As with e, (13), k 12 in (33) may be neglected ab 


ove 


km 


L. 

E,s (6), (7), (9>-(17). and (22)-(32) together with the reaction 
, s in Table 1 comprise the photochemical model to be used in 

Lions 3. A, and 5. 
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3. ENERGY RELATION 
The energy relation is 


dT 
f3C v dt 


= S + pV*^ 


( 34 ) 


where p = pressure, T = temperature, p = density, = heat capacity 


at constant volume. Here 
S 2 E + C 


(35) 


where E represents energy absorbed and C is the heat exchange due to 
infrared radiation. In the atmosphere the fractional pressure change 
is small for velocities less than that of sound (Jefferys 1930). This 


permits ea (34) to be approximated by 


— + W -E- = -J— (E+ Q (36) 

dt c p f c p 

where c « heat capacity at constant pressure, g ** gravitational con- 
P 

stant and w is the vertical velocity. (See Ogura and Phillips, 1962, 
for an analysis of the validity of (31); holton, 1971, shows that (31) 
is rigorously correct in a coordinate system where log p replaces z.) 

As discussed in the introduction the major source of cooling in 
the region 22-61 kra is the 15h band of carbon dioxide. This cooling 
proportional to the local blackuody function (Rodgers and Ralshaw, 
1966) which, for the range of temperatures found in the stratosphere. 
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can be approximated as a linear function of temperature. Thus, the 
cooling or thermal emission, is 


-i— C ~ -a(T - 180°) ... . ( 37) 

. PC P 

where a is a function of height. We use two different distributions of a: 
(1) a = (17 days) - * (Lindzen and Goody 1965) and (2) a = (10 days) , 
z > 35 km. a = (10 days)' 1 (1 + a < 35 km, (Dickinson 1968). These 

two distributions', referred to as slow and fast infrared cooling 
respectively, offer a range of cooling rates. They show the effect tnat 
changing the infrared cooling rate has on the stratospheric temperature 
and ozone distributions. More recently Dickinson (19 72, personal com- 
munication) has suggested that 'a 1 continues to increase above 35 km. 

The energy source may be written as 

£ ■ f * • «t (38> 

where r « is the number of reactions per sec per cm for reaction JL 
and e 4 is the energy gain per reaction. • The £j , consists of dissocia- 
tion energy, D, , and solar energy absorbed, . With these 

definitions eq (33) becomes 

£ =2 ' i +D i >• <39> 

i 

Solar energy absorbed by O y O^, and NO z contributes to the 


source terms: 
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1*L H = hv 0 2 ' 1 2 n 0 2 + h,s 0 3 (q 3a + q 3 b)no 3 


+ hv N0 2 q N0 2 n N0 2 


+ hv H 2 O q H 2 O n H 2 Q 


•where 


q.(z) = Ja.l Q exp( -£ a.x.)d v , j = C> 2# Oj, HO^ NC> 2 


and 


-1 


v = q \ a.vl exp( - 4- a x )dy. 
j J J ov J J J 


(40) 


(41) 


(42) 


2> i. Vi 

Here a.(v) is the absorption cross section in can of the j constituent, 
3 

I is the photon flux incident at the top of the atmosphere, x^. is the 
°v 

total number of molecules per cm 2 of the j constituent between the 

top of the atmosphere and the altitude z, h is the Planck constant and 

.th 

v is the mean frequency of solar photons absorbed. by the j con- 

j 


stituent. Representative variations of (q + q 3b ) n 0 ’ q O-, n O ’ q H O n H~ 

3 Z c 


and a n with height are shown in Figure 1. It is apparent that 

h no 2 no 2 

(%a + l3b- ,B 0 3 y) q 0 2 n O 2 and (q 3_a + V"o 3 > 5 q H 2 0 n H .,0 at al ‘ 
heights. Since hv Q /hv^Q and hv^/hv^ vary between 1 and 1/3, 


the term hv Q (<l 3a + <J 3b > n 0 is much larger than hv H 2 0 q H 2 0 n H.,0 ° r 

h" 0 n at all heights. The latter two terms may be neglected 
°2 °2 °2 . ... 
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in eq (40). Above 30 km h v (q 3a q 3b )n Q >> hv NO q N0"\'O ? ’ 

3 3 c. c. o 

30 km, the ratio (q 3a + ^NO n NO_ ^ 8, provided the mixing 

3 2 

i 

- 8 

ratio of NO_ does not exceed 3*10 (the value used in calculating 

Ct 


Below 


£ ‘M0 2 n N° 2 


for Figure 1). Since Hv- /^ v TvrO 

' 3 2 


0. 7 below 30 1cm, the 

n at all heights and the latter 
N °2 

term is neglected in eq (A0). Then eq (40) becomes 


term hv 3 <q 3a + > ^o^NO^NO., 


r £ E i~ hv O <q 3a' i ' q 3b )l V % 


(43) 


The second part of the energy source term in eq (39) can be written as 




= D(0_)n 


dn /n 
°3 m 


3 m dt 


i D '°2> n m !| S n O /n m + * n (/ 3 V 


+ D(HO_)n 


d "H0 2 / V . , 


2 m dt 


+ 2 D(OH)n m< dt 


, , dn / n 

dn /n HO m 

OH m 2 . 


dt 


dn^ /n 

1 d , OH m 

^>D(H,0)n (~7—n In + 


2"* 2 m dt HO’ m 


dt 


d “ N °2 

) + D(N0 2 )n m -(- r - ) 

m 


(44) 


tVi 

where D(j) is the dissociation energy of the j constituent. We have 
evaluated (44) in detail and have found that 'its contribution is always 
small. We will mercifully spare the reader from the specifics, (N. B. 
(44) is identically zero for equilibrium.) 
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Thus eq (31) becomes 


— — + w -2— = n4> - a (T - 180° j 1 (45) 

dz c 

P . 


where 


" 5 ln, 0 3 (t l3a + ’3b )n n /pc p 


(46) 
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4. RADIATIVE-PHOTOCHEMICAL EQUILIBRIUM 
The photodissociation rates and energy absorption all vary with 
the solar zenith angle and thus local changes in temperature and compo- 
sition have a time scale of a day. The photochemical relaxation time 
scales are less than a day above 35 km for all constituents so calcula- 
tions made with the local time changes in composition neglected should 
still give realistic distributions of constituents. However, the 
thermal relaxation time scale is a minimum of two to three days (see 
section 5). An estimate of the daily variation in temperature due to 
changing zenith angle of the sun can be made. The maximum energy 
absorption per day is at 50 km so the maximum daily .variation in temper- 
ature would be expected at this height. At 50 km the energy absorption 

-3 -3 -1 

rate at noon is about 3 • 10 ergs cm sec . Roughly one half of 
this energy goes to maintain the mean temperature. The noontime flux 
isalso about half the amplitude of the diurnal oscillation in absorp- 
tion. Then, from eq (45), the daily variation in temperature is esti- 
mated 


3(T)d ^ (n<t>)d 

3c ~ pc 

P 


a(T)d 



where the subscript d refers to the daily variation and a = (10 days) -1 , 

—6 *“3 7 l 

P " 10 g cm , and c p ** 10 ergs (g °K) - . Thus Td ~ 2°K which is 

small enough that calculations made with the local temperature change 

neglected should still give realistic estimates of the temperature. 
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In the temperature relation, eq (45), the advcctive terms are 


3T , /3T 

+ w ( — + 

P 


3T , 

u + v 3y ' • 


( 48 ) 


The size of the winds and temperature gradients can be estimated from 
Leovy (1964a and 1964b) , who calculated thermally driven circulations 
in the stratosphere and mesosphere. Between 35 and 61 km the advective 
terms are at least an order of magnitude less than the remaining terms. 

3T . o . n ~ 3 ov ' < ^- o v ^ .6 m sec , 

For example, at 50 km & 3 • 10 K km > dz ~ ^ 


w 0 $ .08 cm 


sec" 1 , and = 10°K Ion" 1 so the vertical advective term 
-6 ... -1 


is 8 ♦ 10~ 6 °K sec" 1 and the horizontal advective term is 

n* 2 • 10~^ °K sec" 1 . Since the cooling terra, a(T - 180°), is 
n, 10 -A sec -1 , the advective terms can be neglected in the temperature 
equation. Similarly, the advective terms in the continuity equations 
can be neglected. For example, at 50 km the vertical scalejor the 

ozone mixing ratio is about 24 km, so the time scale for change due to 

advection is ^ 3 * 10 7 sec. At 50 km the ozone relaxation time is 

5.4 • 10 3 sec so the advective terms can be neglected. Below 30 km, 
however, the ozone relaxation time scale increases rapidly to 
3 • 10 7 sec while the advective time scale decreases. Therefore, below 
30 km the advective terms are important in determining the ozone distri- 
bution in the atmosphere. 

The possibility also exists that the tides excited by solar heating 
could produce mean fluxes of importance to the mean thermodynamic budget. 
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Careful checks indicate that this is not the case below 90 kin. Thus 

radiative-photochemical equilibrium seems appropriate from 35-61 km. 

* 

For radiative -photochemical equilibrium eqs (22) , (23) , and (45) 
become 

7 

0 = C - A<p - B<p\i - FXjp (49) 

2 

0 = D<p + G - Eiji (50) 

0 = -a(T - 180°) + ncp. . -(51) 

The distribution of <p, an d T with height for equilibrium is found by 
solving eqs (49), (50), and (51)- simultaneously. The solution depends 
on the mixing ratios of water vapor and the nitrogen oxides and on the 
values of the rates used for the reactions listed in Table 1. 

As indicated in the preceding paragraphs these three equations are 
valid for the height range 35-61 km. Below this height the photo- 
chemical and thermal relaxation time scales become long enough that the 
neglected terms, particularly the advective terms, should be retained 
to find temperature and ozone distributions which agree with those 
observed. For convenience, eqs (49), (50) and (51) are used for the 
entire height range, 22-61 km, with the stipulation that the resulting 
radiative-photochemical equilibrium distributions are not expected to 
agree with those observed below 35 km. 
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' „ and 

An equilibrium solution also depends on q, q 3a . 1 3 b- %.,C 

which are calculated at each height using -V. (41), (42), and 

The spectral region, 1755-3950 X, is divided into thirt y -eight 

, 7 r n 7000 X into four intervals, 
requency intervals and the region. 4750-7000 A. ^ 

f n O HO, and NO ? is 

Below 61 km nearly all photodissociation of v y z 
du e to absorption of photons in these two spectral regmas. 

.... .» •— "" — *" “ 

,» * - «» — " *“ 

b£ low 3100 X. It has been suggested (Hinteregger 1970) that the 

. -i (1966) for the spectral region 

photon fluxes found in Brinkmann et al. (I960) 

below 1800 X may be too large by as much as a factor of three. 
Below 61 km the absorption of photons with energies corresponding 
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to wavelengths less than 1800 K makes only a small contribution to 
dissociation rates of O^, O^. H^O, anc ^ Therefore the question 

5 G 

of the magnitude of photon fluxes for the spectral region below 1800 A 
is not important for our calculations. 

.th 

Actual cross sections, a., for absorption by the j constituent 

vary with frequency. For each frequency interval the average cross 

section, o’., was used. Cross sections for O are from Griggs (1968) 

3 3 

and Dm and Tanaka ( 19 53). Those for H^O are from Watanabe and 
Zelikoff (1953) and Thompson et al. (1963). Those for NO^ are from 
Hall and Blacet (1952) and Nakayama et al. '(1959). It is more difficult 
to determine the appropriate cross section for in a given frequency 
interval due to the presence of the Schumann-Runge absorption band 
which has numerous narrow rotation lines. One way of treating the 
absorption in the Schumann-Runge .band is to consider the absorption 
cross section for O- as a function of both frequency and path length 

M 

(Hudson et al. 1969). Then it is possible to divide the atmosphere in 
layers vertically and calculate the absorptiqn coefficient of a given 
frequency interval for each layer (Brinkmann 1969). .The atmosphere 
between 61 and 22 km is divided into twenty-four layers. 

The water vapor mixing ratio is not well determined in the 
stratosphere. Some observations (Mastenbrook 1968, Scholz 1970) 
indicate a value -of 2*10 ^ - 3*10 ^ at 29 km and 2. 10 ^ - 6*10 ^ at 
50 km. The value, of the mixing ratio used at a given height is 
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important at lower heights only if water vapor is a significant absorber 
of solar photons or if the hydrogen reactions strongly affect the ozone 
density since ozone is a strong absorber of solar energy. In the height, 
range 22—61 km the attenuation of incoming solar radiation due to 
absorption by water vapor is negligible. It is found that hydrogen 
reactions are dominant in determining the ozone density above 40 km fc* 
a water vapor mixing ratio of 5 * 10 The height above \diich the 
hydrogen reactions are dominant increases with decreasing water vapor 
mixing ratio. Since attenuation of solar energy by ozone is small 
above 50 km, it is only at heights below this that the vertical distri- 
bution of water vapor above is important. The effect of varying the 
mixing ratio of water vapor with height is not investigated here. For 
convenience the mixing ratio for water vapor is assumed to be constant 
with height and is varied from 0, 10 10 10 ,5-10 to 10 for 

supcessive cases. • By this variation we show the importance of water 
vapor in determining the vertical distributions of ozone and hence . 
temperature. 

As indicated in section 2, the vertical distribution for the mixing 
ratio of the nitrogen oxides, NO and NO^,, is not well known. However, 
if downward diffusion of NO from above 70 km occurs, the mixing ratio 
of NO must increase with height. The mixing ratio of nitrogen oxides, 
K, is defined in eq (21), and may be written, by use of eq (14), as 
K = (1 + R 3 )X. 


( 52 ) 
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K is given two different distributions: 


(1) K = 3*10 -8 exp(z - 35) /I-I 

= 3 * 10~ 8 


z ) 3 5 km 
z ( 35 km 


and 


(2) K = 3*10~^ exp(z - 35)/H z ) 35 -cm 
■*-9 z ( 35 km 


(53a) 


(53 b) 


= 3*10 

where H = 9. 25 km. These two distributions are within the range of 
observations discussed in section 2. They should be helpful in 
determining the importance of the nitrogen reactions to the ozone a„d 
, temperature distributions. The terms, high nitrogen and Xo» nitrogen, 

refer to distributions (1) and (2), respectively. 

Reaction rates (5)-(26) are uncertain. Three different sets of 

rates are used for reactions (5) and ( 6 ) (the pure oxygen reactions) 
and these three sets are listed in Table 2. The rates to use for 
reactions ( 7 )-< 26 ) are individually increased or decreased by a factor 
of 10 and the equilibrium distributions of constituents and temperatures 
are recalculated. It should be noted that the rates for reactions 
■ involving IWOj are not known. It is assumed that these rates are the 
' same as those for similar reactions involving H.NO 3 . This assumption is 
admittedly a guess which permits an estimate of the number density or 
HN0 2 . From eq (16) it is found that the number density of K.N0 2> is 
sufficiently small that the equilibrium estimate of UNO., by eq (1C-) is 
adequate for the current work. Even if the actual rates for 
reactions ( 20 ) and ( 22 ) are such as to increase the number density of 
RMO by two orders of magnitude from those estimated in eq (16) using 
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the rates in Table 1, the effect of HNO^ in the current work could 
still be treated by use of eq (16). Since reactions (20) and (22) do 
not appear in eqs (49), (50), and (51), the rates used for these reac- 
tions will not invalidate ' the calculated temperatures or concentrations 
unless the actual rates are such that the number density of is 

more than two orders of magnitude larger than estimated. If the actual 
rates are such that HNO^ number density is more than two orders of 
magnitude larger than estimated here, the validity of eq (16) would 
have to be reconsidered. 

In solving eqs (49), (50), and (51) consideration must be. given to 
the fact that the thermal emission, as represented in eq (51) by 
- a(T - 180°), occurs twenty-four hours a day while the absorption of 
solar energy and photochemistry as given by the other terms in the 
three equations, occur only during the daylight hours. Therefore, in 
eqs, (49), (50), and (51) the terms involving either absorption of solar 
energy or photochemistry are calculated for their noontime values at 
the solar equinox and then these values are divided by two. Above 
50 km where absorption does not vary much over the daylight hours, this 
approximation is fair. However, a.t lower, altitudes it obviously over- 
estimates absorption. That this is the case is shown in time dependent 
calculations which we will mention later; the average temperatures 
obtained from time dependent calculations are generally 10-15°C less 
than the equilibria calculated here. The thermal emission term is used 
at its full value because it operates the entire day. Since the rate of 
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absorption of solar energy at a given height is also a function of latitude, 
equilibrium calculations arc made, at two latitudes: 0° and 45°. 

In Figures 2 and 3 (0° latitude) and Figure’s 4 and 5 (.45° latitude) 
are shown the radiative-photochemical distributions of temperature and 
ozone density for six cases: (1) O (hydrogen and nitrogen reactions 

neglected), (2) H + O (nitrogen reactions neglected), (3) N + O.(low) 
(hydrogen reactions neglected), (4) N + O (high) (hydrogen reactions 
neglected), (5) N + H + O (low), and (6) N + H + O (high). For all six 
cases the rates in set 1 (Table 2) are used for reactions (5) and (6), 
and the fast vertical distribution for a. When the hydrogen reactions 
are included, we use the 5*10 ^ ratio for water vapor. 

Differences among the six different cases are most striking in the 
temperature profile (Figures 2 and 4). The O case shows a linear 
.increase in temperature up to 45 km and then a very small decrease in 
temperature above. The inclusion of the nitrogen reactions (low mix- 
ing ratio for nitrogen oxides) decreases the equilibrium temperature 
and ozone density below 30 km and has little effect on the distributions 
above 30 km. The effect of the nitrogen r eactions is more pronounced 
when the high mixing ratio of nitrogen oxides is used. The temperature 
below 50 km is considerably reduced (<^-'20°K) from the pure O case 
and the ozone density below 40 km is also considerably reduced (^> 40 %), 



31 


Inclusion of the hydrogen reactions has little effect below 40 km 
but above this height strongly reduces the ozone density and tempera- 
ture from the values found for the pure O case. When the hydrogen 
reactions are included, they are dominant in determining the tempera- 
ture and ozone density above 45 km even when the nitrogen reactions 
are also included with either the low or high mixing ratio for nitrogen 
oxides. Most significantly there appears to be no way of producing a 
temperature decrease of the observed magnitude above the stratopause 
without including hydrogen reactions! 

Water vapor must be dissociated before the resultant constituents 
H, OH, and HO^, can act by means of the hydrogen reactions to reduce 
ozone density and thus temperature. It might be expected, since the 
mixing ratio of water vapor is assumed constant with height, that the 
hydrogen reactions would be important at all heights, not just above 
40 km. However, the dissociation of water vapor is strongly height- 
dependent. The photodissociation of water vapor below 61 km is due 
to absorption of solar radiation with energies corresponding to wave- 
lengths below 200Q X, which is strongly attenuated below 85 km due to 

• _ g 

absorption by O- and O . As a result, q decreases from 4.3*10 
« 3 Hy O 

-1 - 11-1 

sec at 61 km to 1.1*10 sec at 22.5 km. Water vapor is also dis- 
sociated in reaction (8) and the rate of dissociation depends on the 
density of O*. From eq (6) the density of O* is proportional to q 
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which is the absorption rate of radiation by ozone below 3100 K. The 
radiation below 3100 R is significantly attenuated below 45 km and q^^ 
decreases from 8* 10 sec at 45 km to 1. 1* 10 sec at 22. 5 km. 

It is apparent that the dissociation of water vapor increases with height 
and, as a result, the hydrogen reactions are more important to equili- 
brium distribution of ozone density and temperature with increasing height. 

In Figure 6 equilibrium temperature profiles are shown for the 
three different sets of rates for reactions (5) and (6) and for the twc 
cooling rate coefficients, fast and slow. The water vapor mixing racio 
used is 5*10 -6 wliile the low mixing ratio for nitrogen oxides is usee. 

The slow cooling rate coefficient (a ^ «* 17 days) is 70% longer than 
the fast one (a ^ ** 10 days) above 35 km. From eq (50) it might be 
expected that a 70% decrease i n 'a* would lead to a 70% increase in 

(T-180°) rather than the 25% increase (at 45 km) indicated in Figure 6. 

. % . 

Such is not the case, since the ozone density itself is strongly temper- 
ature dependent through reactions (5) and (6) and increasing the temper- 
ature decreases the ozone density and thus the absorption of energy. 
Therefore the equilibrium temperature is buffered against change. An 
alteration in the cooling is nearly compensated for by a corresponding 
alteration in the heating so the fractional change in the equilibrium 
temperature is actually much smaller than the fractional changes in 
either the cooling or heating. 
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Coupling of temperature and ozone density is obvious for an O 
atmosphere (neglecting the hydrogen and nitrogen reactions) since ( 49 ) 
becomes 

0 = C-A<p 2 (54) 

and, from eqs (26) and ( 9 ) 



in. which the ratio k^/k,_ is strongly temperature dependent. It should 
be noted that for each set of rates for k c and lc. in Table 2, the ratio 
k^/k,. is strongly temperature dependent. In fact, the temperature 
dependence of k^/k^. for set 1 is exp (-3395/T) which is very close to 
that of exp (-3295/T) found for set 2. The temperature dependence of 
set 3 is somewhat less and the effect of this will be discussed in 
section 5. When hydrogen reactions are dominant in determining the -< 
ozone density (as appears to be the case above 45 km for a water vapor 
•mixing ratio of 5*10 it' might appear that the temperature dependence 
of the ozone density would be greatly reduced, since the hydrogen reac- 
tions do not have rates that are strongly temperature dependent. Above 
45 km eq (49) can be approximated, by use of eqs (27) and (9), as 


0 ** C - B# 


(55) 



34 


where (above 45 km) , 


p ~ Z(q 3a + %b )k 9 


(56) 


It is true that is not temperature dependent but is strongly tem- 
perature dependent. Therefore, the presence of hydrogen reactions 
does not eliminate the buffering of the temperature profile. This 
result should be expected since the hydrogen reactions only destroy 


ozone while the production of ozone by reaction (5) is temperature 
dependent no matter what destruction processes are dominant. Eqs 
(55) and (56) indicate that the ozone density and, thus, temperature 
profile above 45 km are dependent on the rate used for reaction (5). 
The temperature profiles in Figure 6, where three different values 
are used for the rate of reaction (5), show this dependence. 

The mixing ratio of OH may be estimated from eq (50) as 




1/2 


(57) 


The use of eqs (30), (31), and (32) in (57) gives 


4 » = 


^ k 8 q 3b^ /k 7 + q H 2 o'| 

(," " (k 13 + j 



( 58 ) 
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Thus the. mixing ratio of OH is proportional to the square root of the 

water vapor mixing ratio. From eq (49) it is apparent that the impor- 

> 

tance of the hydrogen reactions to the ozone mixing ratio is directly 
proportional to B<pip. Thus, increasing the water vapor mixing ratio 
both increases the importance of the hydrogen reactions at a given 
height and decreases the height at which the hydrogen reactions become 
dominant in determining the ozone mixing ratio. 

Above 45 km eq (13) may be. used for ; i.e., 




n 0 k 2 7 

± \ 

n n n k ' 
C>2 ra 10 



' (59) 


and (58) becomes 


* 


(k 8q3 b<f/k 7 + „>k u n 1/2 „ 



k Q k 1A n n, + —X ) ' 

y 14 m n n k ' 

0 2 m 10 


9 14 in 


Also,, above 45 km, we may use an approximate form of (33); i.e.. 


^ 2R.k Q n 
19m 


Eqs (59), (60), and (61) are not used in our radiative photochem- 
ical equilibrium calculations but they are useful in estimating the 
effect of varying tiie rates for the hydrogen reactions. 
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It is apparent from eq (60) that changes in k g , k^, kg, k^, '^27’ 

and k^ will be ^indistinguishable in their effects from changes in the 

> 

assumed water vapor mixing ratio. The effect of varying the water 
vapor mixing ratio has been discussed previously. We shall return to 
this matter later in evaluating our quantitative results. 

The effect of the nitrogen reactions on the ozone mixing ratio is 
indicated in eq (49) by - F<j)X. Use of eqs (28), (14), and (21) plus 
the fact that R.k,, > k at all heights yields 



and 


F(aX ^ 2R.k,.n Krt below 35 km (63) 

r 1 16 m v 


F (pX ~ Z R 1 k l6 n rn K< p/( 1 + /k 1? ) above 35 km (64) 


It is apparent from eqs (63) and (64) that variations in k^ or k may 
be equivalently expressed as changes in K, the mixing ratio of 
nitrogen oxides. The effect on the ozone mixing ratio due to varying 
the nitrogen oxides mixing ratio has already been indicated. 
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Finally, it should be noted that below 30 km the photochemical 
time scale for ozone is so long (several weeks to years) that dynami- 
cal effects can be significant. It follows that joint radiative-photo- 
chemical equilibrium calculations should not be expecte’d to give 
results in agreement with observed temperature 'and ozone density 
distributions. Above 35 km both the photochemical time scale (( 1 
day) and the radiative-photochemical cooling time scale (about a few 
days) are short enough that dynamical effects are unimportant. There- 
fore equilibrium calculations should give appropriate values of ozone 
density and temperature. 

2 — 6 clearly show that distributions of ozone temperature, 
and density depend on values used for parameters such as reaction 
rates, cooling time scale, and mixing ratios of water vapor and nitrogen 
oxides. These parameters can be adjusted to yield distributions which 
best match the observed distributions of ozone density and temperature. 
Comparison of ozone measurements made by several investigators 
(Johnson et al. 1954, Rawcliffe et al. 1963, Miller and Stewart 1965, 
Weeks and Smith 1968, Ililsenrath et al. 1969) shows a range in number 
density as large as an order of magnitude at some heights in the 35-61 km 
region. Since the temperature in this region is much better determined 
as in the U. S. Standard Atmosphere 1962 or Supplemental Atmospheres 
(Handbook of Geophysics and Space Environments), temperature rather than 
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ozone density is used to determine the match of the equilibrium calcu- 
lations to observations. However, as we shall see, a rather systematic 
difference between our predicted ozone distribution and those observed, 
may indicate an error in our choice of cooling rate coefficient. 

The Supplemental Atmosphere (Tropical) has a stratopause tem- 
perature of 270°K with 243 K at 35 km and 250 K at 61 km. In Figures 
2, 4, and 6 all stratopause temperatures are higher than 270 K. Of 
the three sets of reaction rates for the pure oxygen reactions, set 1 

gives the smallest stratopause temperature. Also, from Figure 6 

-1 ' . “ ' 

the fast cooling rate (a = (10 days) above 35 km) always gives lower 
stratopause temperatures than the slow cooling rate (a = (17 days) ). 

Set 1 for the pure oxygen reaction rates and the fast cooling rate are 
used in the calculations to model the temperature. This choice means 

i <■ 

smaller mixing ratios of water vapor and nitrogen oxides are needed 

to lower the calculated stratopause temperature to that of the Supple- 

' • 

mental Atmosphere (Tropical) than if a slower cooling rate or another 
set of rates for the pure oxygen reactions were used. 

Temperature profiles in Figures 2 and 4 indicate that cases with 
water vapor absent have a nearly constant temperature above the 
stratopause while cases with a water vapor mixing ratio of 5*10 ^ show 

9 

a decrease of 55-75°K in temperature from the stratopause to 61 km. 
Since the Supplemental Atmosphere (Tropical) has a 20°K decrease from 
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the stratopause to 61 km, the water vapor mixing ratio is between zero 
and 5- 10" 6 . A value of 5* 10 _7 gives ^ 23°K decrease between the 
stratopause and 61 km but is less than the 2-6*10 6 value for the water 
vapor mixing ratio that measurements indicate (Mastenbrook 1968, 
Scholz et al. 1970). As eqs (60) and (61). indicate a decrease in the 
dissociation rates of water vapor, a decrease in k R or k^ or an in- 
crease in k would mean the water vapor mixing ratio must be in- 
creased to maintain the same calculated equilibrium temperature pro- 
file. . • - - 

The effect of various nitrogen oxide mixing ratios on the strato- 
pause temperature is shown in Table 3. Profiles 1 and 2 are the low 
and high nitric oxide mixing ratios described previously. Profiles 3, 

-9 

4, and 5 have constant nitrogen oxides mixing ratios of K = 3*10 , 

3*10 -8 , and ITO’ 7 , respectively. Profile 6 has 

K = 3*10”^ 2 (35 km 

■ K = 3* 10~^exp( ~~^~~ ~) 35 km.( z < 45 km 

K = 1. 6* 10” 8 2 ) 45 km. 

Profile 7 is simply Profile 6 with the K values multiplied by 10. Pro- 
file 8 has no nitrogen oxides included. Profiles 2, 5, and 7 yield about 
the same stratopause temperature. Profile 5 with a constant mixing 
ratio is unsuitable since the temperature decreases by 56°K from the 
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stratopause down to 35 km as opposed to a 27°K decrease in the Supple- 
mental Atmosphere (Tropical). Profile 2 has a mixing ratio of nitrogen 
oxides which continues to increase with height above 45 km. Conse- 
quently the temperature decrease from the stratopause to 61 km is 31°K 
for Profile 2 as opposed to 23°I< for Profile 7 which has a constant 
nitrogen oxides mixing ratio above 45 km. Therefore Profile 7 is 
preferred of those tried. The mixing ratio of nitrogen oxides at a 
given height affects the calculations at lower heights only if the 
amount of solar energy reaching the lower heights depends on the nitrogen 
oxides mixing ratio higher up. This dependence occurs if the nitrogen 
oxides are strong absorbers of solar energy or if the nitrogen reactions 
strongly alter the ozone density since ozone is a strong absorber of 
solar energy in the relevant frequency intervals. The former is not 
true and above 50 km, for the mixing ratios of nitrogen oxides used 
here, neither is the latter. Thus, equilibrium distributions found at 
the stratopause depend only on the nitrogen oxides mixing ratio at that 
height. From Table 3, we conclude that, since Profiles 2, 5, and 7 all 
give stratopause temperatures of about the same size, that a mixing 
ratio for nitrogen oxides in the range of 8-16*10 is acceptable for 
the reaction and cooling rates used. However, as will be seen in 
section 6a, our choice of cooling rate may be too low in which case we 
sr<- requiring that ozone density be too low and consequently exaggerating 
the need for ozone los.s processes. A larger value of 'a' will much 
diminish our need for NO . Moreover, we are really specifying FX in 
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(49) rather than the NO mixing ratio, and changes in reaction rates 
would alter our estimate for the latter. 

From eq (64) it is apparent that the mixing ratio of nitrogen 
oxides needed to give a particular temperature distribution is altered 
when the reaction rates, and k^, are altered. 

The treatment of the day-night differences in photochemistry is 
fairly crude. However, time-dependent calculations have been made 
for the daily variations in ozone density and temperature for Profiles 
7 and 8. The temperatures at the stratopause, averaged over 24 hours, 
are listed in Table 3. The stratopause temperature from the time- 
dependent calculations is 15 °- 17°K lower than that found in the joint 
radiative-photochemical equilibrium calculations. For Profile 7 t-|ie 
resulting stratopause temperature of 27 3 K is in good agreement witn 
the 270°K stratopause temperature of the Supplemental Atmosphere 
(Tropical). The results of the time-dependent calculations will be . 
reported separately. 

The 15°K difference between stratopause temperature found for 
radiative-photochemical equilibrium calculations and for the time- 
dependent calculations when the same reaction rates and mixing ratios 
of water vapor and nitrogen oxides are used is due to the crude aver- 
aging technique for day-night differences in the radiative-photochemical 
calculations. As discussed earlier in this section solar energy absorp- 
tion occurs only during daylight hours. The average energy absorption 
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rate was set equal to half of the noontime energy absorption rate. 
was pointed out that this treatment of day-night differences leads to 
an overestimate in the energy absorption rate and, consequently, an 
overestimate in the radiative-photochemical temperature. From eqs (49), 
(50), (51), (55), and (56) it is seen that a stratopause temperature of 
273°K instead of 283°K results when (all' other parameters being the 
same) the noontime energy absorption divided by 2. '4. instead of 2 is 
used for the average energy absorption rate. 

4a, CRITICAL ASSESSMENT. OF EQUILIBRIUM RESULTS 

The above results while interesting are not as important as the 
fact that the system we have developed is simple enough to permit the 
immediate evaluation of the effects of various changes and the identi- 
fication of questionable rate coefficients. In this section we will 
use this simplicity in order to assess and extend the results obtained 
above. 

In the above scheme we have essentially asked what amount of ozone 
is necessary to produce the observed mean temperature between 35 1cm and 
60 km. The amount is less than what one would obtain from simple oxygen 
chemistry. So we further inquired how much H^O and NO^ would be needed 
for the various catalytic destructions of ozone to reduce the value to 
that needed. For the reaction and cooling rate coefficients used we 
obtained answers, to these questions, and these answers seem questionable 
on at least three counts: 

i) The density of ozone calculated between 40 and 50 km is about 30% 
less than the few existing measurements indicate; 
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ii) The suggested mixing ratio for NO^ at 50 km, 1.6*10 , is higher 

than has hitherto been' suggested ; and 

iii) The mixing ratio for 1^0 suggested, 5*10 , is less than currently 

believed. 

For argument's sake, let us assume the above discrepancies are cite 
real, and attempt to track down their origins. Item (i) is entirely 
determined by eq (?» ) ; i.e., it is independent of any chemical assump- 
tions, and depends only on our choice of solar intensity and absorption 
cross-sections in ozone's absorption bands and on our choice of a , 
the cooling rate coefficient. Considerable uncertainty exists con- 
cerning the value of 'a' at 50 km. If it were 30% larger, 301 more 
ozone would be needed to produce the observed T. Dickinson (1972, 
private communication) has informed us that his current calculations of 
15y C0 2 cooling suggest that 'a* indeed continues to increase above 
35 km to values on the order of 30% higher than 1/10 days above '-*0 km. 

With 30% more ozone above 40 km, the loss terms due to hydrogen 
and nitrogen reactions must be reduced. Assuming our reaction rates 
are correct, this can be achieved by reducing our mixing ratios for NO^ 
and H^O. We rule out the latter on two counts: 

a) Our value for n„ n /n already seems low (viz item 3), and, more 

H 2 0 m 

important, 

b) The value we have chosen for n un /n determines not only the value 

' El 

of <p (and hence temperature) at 50 km but also the gradient of temper- 
ature above 50 km where HO has a negligible effect. 
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We are, therefore, forced to reduce our mixing ratio for NO^. To allow 

■’ " -7 

for 30% more ozone, this mixing ratio must be decreased from 1.6*10 
to about 2.4*10~ 8 . The latter value seems more plausible (viz item (ii)) 
and in turn supports the choice of the larger value for 'a'. Thus, a 
plausible choice for the NO mixing ratio between 35 and 50 km appears 
to be 2-3*10 8 . This is considerably more (by about a factor of 20) 
than what is expected will be produced by full scale SST operations 
(Johnston, 1972). Thus, while SST effects should be discernible, they 
should not be catastrophic. In particular, a change of 5/. in NO^ 
should not produce a similar change in <{> since nitrogen reactions .re 
not the only loss processes for ozone and also because of the thermal 
buffering already described. Turning, finally, to item (iii) we must 
note that the quantity we have chosen is not really n fi Q /n m but B\j/ in 


eq (49). From eqs (60) and (61) we have that 




2(, >3a + q 3b )k 9 
\ k 5 


1/2 1/2 

(k 8 q 3b 4,/k 7 + q II o 0 )k ll 1 n H 0 

- -r— P : '- 

n 0 k 27 -'■■■ ' m 

k n k . , n ( 1 + — r — ) J 

9 !4 m n 0 2 m^lO 


above 45 km. 

Clearly, various changes in reaction rate coefficients could allow 


n 


H 2° 


2 Vciluc for ( — * *) without clisnging RcCcnt messurernents 

in 

brought to our attention by S. Wofsy suggest our rate for k g should be 
increased from S-IO” 11 to 2*10 _1 ° while k^ should be changed from 
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10 -11 to 2*10~ 10 (see Wofsy et al. 1972 for references). With these 

n H 2 ° _ 6 

new rates r.ust be increased to about 3-10 to maintain a given 

n 

m 

value for By. Discrepancy (iii) argues in favor of the new rates. 
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5. RELAXATION TIME SCALES 

Suppose we consider small perturbations from the radiative- 
photochemical equilibrium state calculated in section 4. Then lineari- 
zation of eqs (22), (23), and (45) yields the equations which can be 
solved for the time scales for return to radiative-photochemical equi- 
librium. The equations are, with winds neglected, 

-r^' = -ft to' OJi' T', (65) 

( 66 ) 

. (67) 



where, using primes for perturbation quantities and overbars for 

9 

equilibrium value, 

<p = <p + <p' , 

= ijj + ip* , 

T = T + T\ 
n i t 1 

m - — (valid for small pressure perturbations) 

' A* ~ 

- T 

n 

m 
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*8 = [2A(p + Bi|i(l + 6 X ) + FX(1 + Gj) J /(I + R x ) 
C = Bp/(1 + Rj), 

^ “ ''2 *2 

= [D - gE?/^]/(l +.R 2 X 
£ = 2 E ^/{1 + R 2 ), 


( 68 ) 

(69) 


& =^[<*A^ + (p - 6 2 g 2 )B4/ + (f + e 2 )FX]/(l + R x ), (70) 
~~ 2 • (71) 

(72) 


and 


J = [rKif> 2 /T]/(l + R 2 ) . 


Further definitions are 


<*= 1 + e 6 - e gl 
Y 4 = Yi + Y 2 * 


(73) 


Y 1 


k l6 k 15 


' 17 * 


R 3 Xk 15 


x/p^ + 


R 


1 * 


k l5 Xq NO, 


V 2 = ~~ — 2— , _ /[k ll + 


R 3 Xk 15 




R i* n m k 17 


R 


1 * 


k i2 + Rjkg^lZ " (2 ‘ e 5 ) - ] ‘ Y l (e l5 + e i6 ' S 17 ) 


“ 7 2 [ e 15 + 1 - e 17 - ( 2 “ e 5 )]» 
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VSno/P + R 3>(<P" m V^ 

Ct 


R 1 1c 16 


e 2 =-(l + R 3 )k i7 


[ e i6 + {Z ~ e 5 ) ' e l7 ] ‘ e l (1 " e l7 ) ’ 


g “ Y* z < 40 km 


2 ^9 ' 

g - -A/ (1 + A) • j— z > 40 km 


A = <J> 


27 


n o 2 k io 


y 3 » Y 4 (l + e^/R^ 


r - w - c 2 y 4 - 1 


z < 40 km 


’ A . -JL ( e + i) + — 1 jL- ( e - 2 + e ) 

(1 + A) 2 k i£*i 27 Vii^ 12 3 


z > 40 km 


£ = - 1 + [2Rjkj 6 (e J6 + 2 - e 5 ) + + kjg] 


p = - 1 + [{k 9 + R 2 k n )(2 - c 5 )Rj + RjRgRjja; 


+ k e ]n /B 
12 12 J m 
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6 , = 

1 '3 1 Z 11 m 


6 2 = 6 lV Y 3 


- dk. 

e i = k7 ( dF ) f ' i = 5 - 18 >- 

and, from the linearization of eq (24), 


X' cp 1 , T 1 

~ = e l- + c 2~ 
X <p T 


Hie complexity of our coefficients in (65) and (66) arises because, 
some of the reaction rates, k^, and the ratios R^, R ? , and are tem- 
perature dependent and the ratios and R^ are also dependent on the 
ozone mixing ratio. There are additional terms in eqs (65) and (66) 
which have been neglected since they are much smaller than those 
retained. If there is no coupling, then the time scale for return to 
equilibrium is ^ for the ozone mixing ratio, £ ^ for the hydroxyl 

mixing ratio, and a ^ for the temperature. The terms Cy' , $1 T ' , <j) ' , 

■^T' ,. and H ^ * > represent the coupling between the photochemical and 


radiative processes. 



50 


Eqs (65), (66), and (67) may be combined to: 


' [ 

where 

M = a + t3 + £ 

N=l3£ + al3 + £*-{- )1 n+ C D 

S = £a$ + + Ca J Cfi • 

Initial conditions are needed to solve eq (74) and we consider three 
cases: . - - 

(1) <p'(t = 0) = <p(0), 

4i'(t = o) = o„ 

T'(t = 0) = 0, 

' (2) <p l (t = 0) = 0, . 

^>(t = 0) = 4(0), 

T'(t = 0) = 0, 
and (3) <p'(t = 0) = 0, 

4»'(t = 0) = 0, 

T'(t = 0) = T(0). 

The solution for <p' for case 1 is 


— o + N-^+ S])J\= 0 

bt fit * IT' 


( 74 ) 


<p' = cp(0) K u e 




+ K 21 G 


°2 t 


+ K 3i e 


a 3 t 
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the solution for for case 2 is 

o,t . o 2 t o 3 t 

= +(0) K 12 e + K 22 e + K 32 e . 

and the solution for T 1 for case 3 is 

0,t • O-t 0,t 

T' = T(0) K^e + K 23 e J + K 33 e 


The coefficients are defined as 

i 

/ 

K, . = ( 0 . + V.(o 2 + a 3 ) + o 2 o 3 )/[(a 1 - o 3 )(a 1 - o 2 )]; 

3 3 3 

K 2‘ = + V '^ a 3 + °1^ + °i°3^L(o 2 - o 3 )(a 2 - o^]; 

3 3 3 

K 3 . = (0^ + V.(c^ + o 2 ) + °i° 2 ^/^( a 3 ” a i^°3 ” a 2^ 
•where j = 1, 2, 3 and 

v L = 13. 

v 2 = X. 
v 3 = <k > 

' Bj = -B 2 - n^3 - CA 
0 2 = £ 2 -^. 

0 3 = a 2 - . . 

The values 0 , 0_, and 0_ are roots of the cubic equation: 

i. Lt -5 

0 3 +/io 2 + No + s = 0 (75) 


which follows from eq (74). 
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For cases where the hydrogen reactions are neglected (74) reauces 


to: 


[ 7 + ( a + $ )~+ (d ^ + = 0 


2 
st . 


St 


(76) 


with the solution for <p' (case 1) as 


<P 


1 — 


= (p(0) 


o 2 +t3 ojt e 2 t 


°2 " °1 


°2, 


(77) 


and the solution for T' (case 3) as 


T' = T(0) 


o 2 + Ot o^t + 0. o^t 
e - e 


A 

Cb 


A 

O-, 


°2 " °1 


(78) 


where 


\ = \ - (a +13) - [(a -^) 2 -.4niJ ] 


1/2 


(79) 


0 2 = "2 ■ ( a ) + ■[( a -*$) - 4^-S] 1/2 . 


(80) 


In Figure 7 time scales a' 1 ,# 1 , T } (= - ( Re^ )’ 1 ), and r 2 (= -(Re^)‘ 

are plotted for the pure O case with set 1 for rates of reactions (5) and 

(6) and the fast profile of a. The possibility of complex a and c 

1 2 

merely means that the coupling results in an oscillating decay of the 



53 


perturbations. The time scale ~13 is several years at 22 km and 
decreases to about 90 minutes at 50 km and then increases above to 
about 1/10 day at 6 1 km. When coupling is included the actual time 
scale for thermal relaxation is T ^ an d for ozone relaxation is iv,. As 
shown in Figure 7 coupling results in a slightly more rapid ozone relax- 
ation below 30 km and a slightly slower ozone relaxation above 30 km. 

Below 30 km coupling causes a slower thermal relaxation. If we 
write eqs (79) and (80) as 


Oj = ” ( a + ) - 6/2 


(81) 


a 7 = - ( a +13 ) - 6/2 
2 


;s2) 


where 6= [( a -"23 ) - 4$ ri] and then substitute eqs (81) and (82) 

into eq (78), the result is 




T' = ^—e' 2 [ ( 13 - a )(e 

. 46 


6 /2t -6/2t 


6 /2t _ 6/2t 

+ 6(e + e )]. 


(S3) 
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Now 5/2 ( ( Cl + 3 )/2 and if wc consider an expansion 
( (6/2) ^ then we may write eq (83) as , 


of 


+ 6/2t 


for time 


i-S4l.lt (g-fl) 

T' £ T(0)e" [1 - 2 ] (84) 

or 

T' ^ T(0)e“ a t . ( 85 ) 

Thus, if we are considering a- wave with period r < <2 ** then the actual 

time scale for thermal relaxation is not 7 -. but (X below 30 km where 
eqs (84) and (85) are good approximations' to eq (83). - 

Above 30 km "ft " 1 is much less than a" 1 , and coupling between 
the photochemical and radiative processes results in more rapid thermal 
relaxation as shown in Figure 7. Lindzen and Goody (1965) also found, 
that thermal relaxation was more rapid above 30 km, though they 


used a different set of rates for reactions (5) and ( 6 ). 

Above 40 km, 4(& p+al3) (< (a + 13 ) and 13 ' ((a" so eq 

(81) can be approximated: 


A & h 

°1 ~ " ■£ . (86) 


For the pure O case eqs ( 68 ) and (70) reduce to 

= 2A<p 

=C4A<pf. 

T 


(87) 

( 88 ) 
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Use of eqs (87) and ( 88 ) in eq ( 86 ) yields 


A c-cA-tco/T C’s'A’.O 

Oj ~ - ZZ _ • 


2 A, 


<P 


2T 


( 89 ) 


Substitution of eq (51) into eq (89) gives 


A ~ Q< a 


a. 


180 

(1 z“). 


( 90 ) 


From Figures 2 and 4 we see that the equilibrium temperature for the 
pure O case is nearly constant above 40 km so, from eq (90) 


~ constant above 40 km 


as»is shown in Figure 7. 

The N + H + O case requires the solution of eq (75). In Figure 8 


as well as 7 -^, t 2 > anc ^ 


O' 1 , "8 and g' 1 

T, (= - [Rec^]" 1 ). where is the hydroxyl relaxation time scale when 


are plotted the values 
,-l 

3 v- - l — v 3 

coupling is included. Set 1 is used for rates for reactions (5) and ( 6 ) 


and the fast profile is used for <2 . The water vapor mixing ratio is 
5. 10~k w hile the low mixing ratio profile for the nitrogen oxides is used. 
In Figure 8 we see that below 40 km 7 - 3 ~ £“ 1 and ^ and r 2 are approxi- 
mately the same as for the pure O case. This result is reasonable since 
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-8 - 1 


we previously found that inclusion of the hydrogen reactions has little 
effect on the equilibrium distribution of temperature and ozone density 
below 40 km. The only effect of the nitrogen reactions is to shorten 
at the lowest altitudes. 

Above 40 km the hydrogen reactions become increasingly important 
with £ ^ nearly the same size as 13 \ The time scale for relaxation 

of a perturbation in the ozone density or hydroxyl density is modified since 
there is strong coupling between the ozone and hydrogen photochemistry. 
The photochemical acceleration of the thermal relaxation is modified 
but not eliminated. Above 40 1cm eq (86) is applicable and: 

Ojfc -JtWB. (91) 

For the N + H + O case above 40 km eqs (68) and (70) reduce to: 

13 & Bvji, • ' (92) 

ft ~ PB^/T (93) 

so substitution of eqs (92), (93), and (53) into eq (91) yields 


* 18 Cj 

Prrp/T = pa(l - — — ). 


(94) 


Now the equilibrium temperature for the N - H - O case in Figures 2 and 
4 decreases nearly linearly above 42 km and this decrease is reflected 
in the nearly linear increase of with height in Figure 8. 
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Figure 9 shows 7 ^ for six cases: O, H + O, N + O (low), N + O 
(high), N + H O (low), and N + H + O (high). It is apparent that inclu- 
sion of nitrogen reactions for the low mixing ratio of nitrogen oxides 
does not alter 7 from cases with no nitrogen present but use of the hi-ti 
mixing ratio increases 7 ^. The reason for this effect can be seen mo.-/, 
readily by referring to the equilibrium eqs (49), (50), and (51). The 
effect of nitrogen reactions on ozone density is contained in the term 
- F<pX and above 35 km this term can be written, by application of 

eqs (29), (14), and (9), as: - ” ' ' , 

FcpX ^ 2k 1? K<pn m (95) 


The effect of the pure oxygen reactions on ozone density is contained 

2 

in the term, - A<p and from eqs (26) and (9), we have 


2 2 (t l3a + q 3b> 


A <p = 


n. 


(k 6 /k 5 )<p 


(96)' 


Comparing the terms (95) and (96), we see that temperature dependence 

of ozone density due to nitrogen reactions is by k^ and the temperature 

1/2 

dependence of ozone density due the oxygen reactions is by • 

1/2 

The temperature dependence of k I s touch less than that of (I^/^cp * 

When the high mixing ratio profile of nitrogen oxides is used, F<pX 
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becomes more important than Ajp in determining the ozone density 
between 35 and 45 km. Since the temperature dependence of F<pX is less 

"I 

than that of A<p , it follows that 7 is increased at heights where nitrogen 
reactions are more important than pure O reactions in determining ozone 
density. 

Figure 10 shows for the N 4 H f O (low) case when three dif- 
ferent sets of rates for reactions (5) and ( 6 ) and two values for CL, 
fast and slow are used. Above 35 km eq (94) still holds, 

o Z PA(1 - i5 ±-)'. .... (v?) 

T 


Since the equilibrium temperature varies with different sets of rates 
for reactions (5) and ( 6 ) (Figure 7), so does 

As discussed in section 4, the effect of varying rates for reactions 
(7)- (18) can be shown by varying the water vapor mixing ratio or the 
mixing ratio of nitrogen oxides. Resulting values of 7 are within the 
range of those shown in Figures 9 and 10. 

In section 4 we found that the equilibrium temperature calculated 
had the best match with the observed temperature for the case with 


set 1 for the rates for reactions (5) and ( 6 ), the cooling time scale cf 

_7 

O =10 days, a water vapor mixing ratio of 5-10" , and a nitrogen 
oxides mixing ratio defined by Profile 7. Using these parameters in 


■the present calculations we find that the minimum value of 


T 1 


is 
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2.6 days at 32 km. Between 35 and 61 km the value of is between 
3.3 and 3.7 da' f s. Since is proportional to 'a', the revisions 
suggested in section 4a will lead to further reductions of T^. 

The results show that coupling between radiative and photochemical 
processes leads to an acceleration of thermal relaxation above 30 km 
for a wide range of reaction rates and other parameters. The minimum 
value of r x is between 2 and 4 days at 31-35 km. When the hydrogen 
reactions are neglected t is nearly constant with height above 40 km 
and when the hydrogen reactions are included with a mixing ratio of 
5*10 _ ^ for water vapor t increases almost linearly from the minimum 
value at 35 km to about 5 days at 61 km. Therefore the Newtonian cool- 
ing rate coefficient to be used in the temperature equation for dynamical 

problems in the stratosphere above 30 km is not 'a* but \ 

seems large enough to dissipate internal Rossby waves. 


Moreover 
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6. CONCLUSIONS 

In this paper we have developed quantitatively justifiable simple 
models for the radiative heating and the complicated ozone photo- 
chemistry of the stratosphere between about 20 and 60 km. These models 
permit the immediate evaluation of the effects of various chemical and 
physical quantities on the temperature and ozone distributions in the 
stratosphere. We show that the region 35-60 kn should (away from 
arctic regions) be in approximately joint radiativc-photocheroical equi- 
librium, and we use our simple models to calculate this equilibrium. 

We show that the interaction of radiation and chemistry significantly 
buffers the temperature and ozone distributions in the stratosphere so 
that rather large changes in cooling rates, reaction rates, etc., lean 
to much more modest changes in ozone and temperature. We also show 
that given an observed distribution of mean temperature (a quantity 
much more extensively measured than ozone) we can (assuming accurate 
knowledge of reaction and infrared cooling rate coefficients) determine 
the amounts of water vapor and nitrogen oxides needed to produce the 
ozone distribution which would in turn lead to the observed temperature 
This procedure is facilitated by the fact (which we demonstrate) that 
hydrogen reactions are significant loss processes for ozone above 40 km 
while nitrogen reactions are important below 50 km. Thus there are 
regions where nitrogen and hydrogen reactions are individually rather 
than jointly important. Unfortunately, reaction and cooling rate coef- 
ficients are not accurately known. However, our simplified equations 
can still be used as diagnostic tools. In .this manner, we show that an 
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internally consistent description of observed distributions of temper- 
ature and ozone in the stratosphere can be obtained with. a time scare 

> 

for infrared cooling varying from ten days at 35 km to 7.7 days above 

~ 8 

40 km, a mixing ratio for nitrogen oxides (NO + NO^) of about 2-3*10 

-~6 

and a mixing ratio for water of about 3*10 . The question of normal 

NO mixing ratio has assumed considerable importance recently because 
x 

it is believed that proposed SST (super sonic transport) operations may 
add about 1*10~ 9 to the NO mixing ratio. The present calculations 
suggest that the SST contribution will be only 3-5% of the ambient 
amount, and because of the thermal buffering of the system the effect 
of this addition on ozone density should be less than 3-5%. 

We also used our simple models to calculate the joint radiative- 
photochemical relaxation of perturbations in temperature, ozone mixing 
ratio and hydroxyl mixing ratio away from their equilibrium values. We 
find that for all photochemical models considered, the photochemistry 
greatly accelerates thermal relaxation as calculated simply on the 
basis of infrared cooling. The time scale for the latter is about 
10 days. The actual time scale for thermal relaxation when photo- 
chemistry (which is temperature dependent) is included proves to be 
only 2-4 days. These shorter scales imply greatly enhanced damping for 
hydrodynamic waves in the stratosphere. 
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Table 1 



Reaction 

Rate 

Source 

1. 

O z + h v —yO + O 

q 2 ( X < 2460 K) 

See text 

2a. 

C> 3 + hv *->0 + 

q (3100 £ < X < 11000 £) 
3s 

See text 

2b. 

+ h v O* + 

q Jb (X < 3100 K) 

See text 

3. 

H^O + hv H + OH 

q H 0 (X <240 ° ^ 

See text 

4 . 

NO^ 4 h v — >NO + O 

' ■ q NO^ < 3975 

See text 

5. 

O + O, + M — >0 + M 

k 5? 

See Table 2 


2 3 

> see Table 2 


6. 

O + 0 3 ~* 20 2 


See Table 2 
1 

7 . 

O* + M O + M 

k ? = 5-10' 11 

Zipf 1969 

8. 

O* + H 2 0 —3^ 2 OH 

k 8 = 5- 10*“ 

Demore and Raper 1967 

9. 

OH + O -> H + 0 2 

kgrrS'lO' 11 

Kaufman 1969 

10. 

H + O z + M -^H0 2 + M 

k io = 4 - 10 ' 32 

Kaufman 1969 

11. 

H0 2 + O -> OH + 0 2 

k^ =5*10 

Nicolet 1970 

12. 

OH + 0 3 — ^HO z + 0 2 

k ? 1. .^10" 12 T 1/Z exp(-1500/T) 

Nicolct 1970 

13. 

oh + oh -m 2 o + O 

k 13 = 5 ' lcrl3Tl/ ^ ex P(- 1000/T) 

Nicolet 1970 


Table 1 (Cont'd. ) 



Reaction 

Rate 


Source 

14. 

OH + HO z ~^r H z O + O z 

k 14 = 

10- U 

Kaufman 1969 

15. 

H0 2 + NO — 5*OH + N0 2 

k 15 

-12 1/2 

3*10 T exp(-l250/T) 

Nicolet 1970 

16. 

NO + O > NO 4 O 

L* L* i 

k l6 = 

3. 2*10‘ 11 exp(-300/T) 

Klein and Herron 1966. 

17. 

NO + 0 3 -^> NO z + 0 2 

k 17 = 

1. 7‘10" 12 exp'(-1310/T) 

Schofield 1967 

18. 

N °2 + °3 ^ N °3 + °2 

k 18 = 

10" n exp(-3500/T) 

Johnston and Yost 1949 

19. 

NO + OH + MH>HNO + M 

k 19 = 

3.3«10" 33 

.Berces and Forgeteg 1970 

20. 

NO + OH + M— *HN0 2 + M 

k 20 = 

k 19 

See text 

21. 

HN0 3 + O NO s + OH 

k 21 = 

1.7* 10' 11 

Berces and Forgeteg 1970 

22. 

HN0 2 + O — P OH + N0 2 

k 22 = 

k 2i . ; . • 

See text 

23. 

NO + NO— > 2NO 

k 23 = 

10' 11 , 

Berces and Forgeteg 1970 

24. 

ho 2 + ho 2 h 2 o 2 + o 2 

k 24 = 

5.10" 12 exp(-1000/T) 

Nicolet 1970 

25. 

H 2 0 2 + h V 2 OH 

X ( 5650 K 

Holt 1948 

26. 

H ? 0 + OH — : >11 O + H0 2 

k 26 

4. 10" 13 T I/2 cxp(-6Q0/T) 

Greiner 1968 

27. 

il + 0 3 • — > OH 4- 0 2 

k 27 = 

1,5-W- 12 I 1/2 

Phillips and Schiff 19^2 
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Table 2 


Set 1: 


k = 1. 2-10" 35 exp(1000/T)~) 


5 , 


.-11 


Schiff 1969 


k 6 = 2. 0* 10 exp(-2395/T)J 

Set 2: 


Benson and Axworthy 1965 


k = 8.0-10 _35 exp(445/T) 
k, = 5. 6*10“ n exp(-2850/T)\ 

6 J 

Set 3: - 

k 5 = 5. 5.10- 34 (300/T) 2 - 6 l Kau , manl967 
k, = 1.4‘10' 12 exp(-1500/T) j 

6 J 



Table 3 




3 4 5 6 7 8 

3*10*^ 3*10 -8 IQ' 7 1.4*10 =8 1.6*10 7 0 

306 301 291 303 288 307 

273 290 

*« 

» 

• ’ I 

i 

1 t 

\ 
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Figure Legends 

Figure 1. The absorption of solar photons by NO^, O^, H^O, and 
as a function of height. 

Figure 2. The radiative- photochemical equilibrium temperature at 
0° latitude as a function of height. The_ six cases are described 
in section 4. 

Figure 3. The radiative-photochemical equilibrium ozone number 
density at 0° latitude as a function of height. 

Figure 4. Same as Figure 2 but at 45 latitude. 

Figure 5. Same as Figure 3 but at 45 latitude. 

Figure 6. The radiative-photochemical equilibrium temperature as a 
function of height for 0° latitude for an N + H + O atmosphere. 


* The three sets of reaction rates are in Table 2 and the fast and 
slow cooling rates are discussed in the text. Our low NO^ values 
were used. 

" 1 

Figure 7. The relaxation time scales in a pure O atmosphere. U 
and a” ^ are the relaxation time scales for O^and temperature, 
respectively, when the coupling between the photochemistry and 
temperature is neglected. t 2 and Ti are the relaxation time scales 
for the an <3 temperature, respectively, when the coupling is 


considered. 
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Figure 8. The relaxation time scales in an N + H + O atmosphere. O , 

C , an d a are the relaxation time scales for O^, OH, and 
temperature, respectively, when the coupling between the 
photochemistry arid temperature is neglected. T,, , 7y and T ^ are 
the relaxation time scales for Oy OH, and temperature, respectively, 
when the coupling is included. 

Figure 9. The thermal relaxation time scale, as a function of 

height for the cases considered in Figures 2 and 3. 

Figure 10. The thermal relaxation time scale, T ^ , as a function of 
0 

height for the cases considered in Figure 6. . 
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